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The double barrier resonant tunneling diode exhibits complex spatio-temporal patterns including 
low-dimensional chaos when operated in an active external circuit. We demonstrate how autosyn- 
chronization by time-delayed feedback control can be used to select and stabilize specific current 
density patterns in a noninvasive way. We compare the efficiency of different control schemes involv- 
ing feedback in either local spatial or global degrees of freedom. The numerically obtained Floquet 
exponents are explained by analytical results from linear stability analysis. 
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I. INTRODUCTION 

Since the groundbreaking work by Ott, Grebogi and 
Yorke Q, chaos control has evolved into a central is- 
sue in nonlinear science 0. While earlier methods 
of chaos control have used a rather complicated cal- 
culation of the control force from the Poincare map, 
recent control schemes based on time-delay autosyn- 
chronization 0, 0] are much simpler to handle and 
have been applied to a number of real world problems 

010110000. 

One intriguing aspect is the possibility of noninvasive 
control. This refers to the stabilization of a target state 
which is not changed by the control term, and the con- 
trol force vanishes once the target state has been reached. 
A natural choice for the target state are unstable peri- 
odic orbits (UPOs), since they are dense in the chaotic 
attractor of the uncontrolled system. 

While earlier work has concentrated on low- 
dimensional dynamic systems described by maps or ordi- 
nary differential equations 0] , the emphasis has recently 
shifted towards stabilization of spatio-temporal patterns. 
It was shown that for a generic nonlinear reaction- 
diffusion model of activator-inhibitor type with one spa- 
tial degree of freedom, different noninvasive time-delayed 
feedback methods can be used to suppress chaotic behav- 
ior 0,0,0,0, and their respective domains of control 
have been compared and interpreted in terms of Floquet 
spectra. 

In this paper we will apply time-delayed feedback con- 
trol schemes to a semiconductor nanostructure which is 
currently of great interest 0|: the double barrier reso- 
nant tunneling diode (DBRT). This device is well known 
as an electronic oscillator, and complex spatio-temporal 
patterns of the current density have been reported in nu- 
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FIG. 1: Schematic energy band structure of the DBRT. The 
symbols are explained in the Appendix. 

merical simulations 0, l2ll|2l. 0, 0] , including chaotic 
spatio-temporal scenarios |2 r l . Here we propose to apply 
time-delay autosynchronization to stabilize those spatio- 
temporal breathing and spiking patterns under a wide 
range of operating conditions, and induce stable periodic 
oscillations. The paper is organized as follows. After 
the introduction of the DBRT model in Sect. [H] we ex- 
amine the dynamical bifurcation scenarios leading to the 
formation of lateral current density patterns of the un- 
controlled system in Sect. lIIII We find parameter regimes 
featuring chaotic breathing as well as chaotic spiking. In 
Sect. lIVI time-delaved feedback control is used for the sta- 
bilization of periodic spiking and breathing oscillations. 
We compare the control performance for different con- 
trol schemes. The details of the model are given in the 
Appendix. 



II. THE MODEL 
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The DBRT is a semiconductor nanostructure which 
consists of one GaAs quantum well sandwiched between 
two AlGaAs barriers along the z-direction (cf. Fig. QJ. 
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FIG. 2: (a) Nonlinear kinetic function f(a, u) for different val- 
ues of voltage u, (b) diffusion coefficient D(a) as a function of 
the electron density a (double- logarithmic scale) , (c) current- 
voltage characteristic for spatially homogeneous states (solid 
line) and filamentary states (dashed line) for L = 30; load 
line for Uo = -84.2895, r = -35 (dotted line), (d) spatial 
profile of a half filament, a(x) (solid line) and j(x) (dotted 
line) (L = 30; u = 265) 



The quantum well defines a two-dimensional electron gas 
in the x — y plane. The spatially homogeneous steady 
states give rise to a Z-shaped current -voltage character- 
istic [26| exhibiting bistability in a range of applied volt- 
ages. The middle branch of the current-voltage char- 
acteristic can be stabilized by applying an appropri- 
ate active circuit with a negative load resistance |27| . 
Complex spatio-temporal patterns can arise if the lat- 
eral re-distribution of electrons within the quantum well 

[SllillilliJlUlHii] is taken into account - In tlic 

following we assume that the extension of the device in 
y-direction is small and charge inhomogeneities can only 
appear in the ^-direction. Using dimensionless variables 
throughout we arrive at the following equations: 



dm 
~dt 



1 



(U -u-r(j))-KF u (t) 



(2) 



Here the uncontrolled model [25| has been extended by 
the control terms KF a and KF U representing the con- 
trol forces with amplitude K. The dynamic variables are 
the inhibitor u[t) and the activator a(x,t). The one- 
dimensional spatial coordinate x corresponds to the di- 
rection transverse to the current flow. In the semicon- 
ductor context u{t) denotes the voltage drop across the 
device and a(x, i) is the electron density in the quantum 
well. The nonlinear, nonmonotonic function /(a, u) de- 
scribes the balance of the incoming and outgoing current 
densities in the quantum well (Fig. 0(a)), and D(a) is an 
effective, electron density dependent transverse diffusion 
coefficient (Fig. 121(b)) [2!j. The local current density in 
the device is j(a, u) = \{f(a, u)+2a), and (j) = j; J jdx 
is associated with the global current. Eq. @ represents 
KirchhofFs law of the circuit in which the device is op- 
erated. The external bias voltage Uo, the load resistance 
r, and the time-scale ratio e are dimensionless external 
parameters. We consider a system of width L with no 
charge transfer through the lateral boundaries (i.e., Neu- 
mann boundary conditions d x a = at x = 0, L). In 
Appendix lAl we give explicit expressions for f(a,u) and 
.D(a) and draw the connection to the microscopic physi- 
cal parameters. 

In the case without control force (K = 0) one sta- 
tionary solution of Q is given by the homogeneous so- 
lution a{x) = a} om with f(a% om ,u) = 0. Up to three 
different solutions a^ om may exist for one fixed value of u 
(Fig-GI a ))- This gives rise to a Z-shaped current -voltage 
characteristic j(u) = a (solid line in Fig. IHc)). 

The spatially homogeneous stationary solutions of the 
coupled equations and @ are given by the intersec- 
tion of the current-voltage characteristic j(u) and the 
load line which is the nullcline (j) = (Uo - u)/r of 
States on the middle branch of the j(u) characteristic 
are unstable in a passive external circuit with effective 
resistance r > 0. By choosing r < 0, which can be re- 
alised by an active circuit, i.e. applying an additional 
control voltage proportional to the device current (j) in 
series with the bias Uo H3, it is in principle possible to 
stabilize the middle branch of the stationary j(u) charac- 
teristic and access it experimentally, but for large enough 
e Hopf bifurcations can occur, leading to uniform limit 
cycle oscillations of the current and voltage. 

It should be noted that, depending upon the intersec- 
tion of the load line with the j(u) characteristics, the 
DBRT represents either a bistable, or an oscillatory, or 
an excitable active medium ,30]- It is remarkable that 
the DBRT can be operated as an excitable system even 
in the spatially homogeneous case. Consider the situa- 
tion in Fig.^fa) where the load line intersects the homo- 
geneous characteristic near the second turning point at 
a low current value (j ps 4.63, and e = 70). Due to the 
large value of e the relaxation in u is much slower than 
the relaxation in a. If noise pushes the system off its 
fixed point in the low current regime, it quickly relaxes 
to a value of a = j in a high current state. Then the 
slow relaxation in u sets in and drags the system along 
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FIG. 3: DBRT in the excitable regime: (a) Spatially homo- 
geneous phase portrait of current j = a vs. voltage it for 
r = -2, U = 248.11, L = 30, e = 70. Solid line: trajectory 
of the system under noise. Dashed lines: Load line (straight) 
and homogeneous current -voltage characteristic (^-shaped). 
Dotted lines: phase flow of the system, (b) Response of a and 
u triggered by noise. 



the homogeneous characteristic towards high u (since we 
are above the load line) until the upper branch of the 
homogeneous characteristic ends. At this point, the sys- 
tem again quickly relaxes in a, this time towards a low 
current state. Since we are now below the load line the 
system finally returns to the fixed point. This excursion 
in phase space results in a spike-like response in the cur- 
rent and voltage signal upon external noise (Fig. Etb)) 
which is characteristic for an excitable medium. 



III. SPATIO-TEMPORAL SCENARIOS 



As the device width L in x-direction increases (the 
device width in y direction is fixed and small compared 
to L), the middle branch of the homogeneous solution 
ghom (Fig.|2Ic)) becomes unstable against inhomogenous 
fluctuations. A straightforward linear stability analysis 
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FIG. 4: Spatio-temporal breathing patterns: electron density 
evolution, phase portrait, and voltage time series for (a) e = 
7.0: periodic breathing, (b) e = 9.1: chaotic breathing; load 
line as in Fig. [5] 



19] shows that this occurs at 



L > TT\ 



D(a\ om ) 



d a f(a%° m ) 



(3) 



Then an additional stationary solution of eq. Q arises 
in the form of a current filament. Fig. |2{d) shows the 
dependence of a and j on the transverse spatial coordi- 
nate x for a half-filament with L — 30 and u = 265. 
The dashed line in Fig. |5fc) depicts the corresponding 
current-voltage characteristic (j)(u). 

The filament is stable for r < at small e, but becomes 
unstable by a Hopf bifurcation for large enough e. In 
this case complex spatio-temporal breathing (Fig.^J and 
spiking patterns (Fig.[SJ are expected if the homogeneous 
fixed point is still stable with respect to Hopf bifurcations 

m- 

Choosing the load line as in Fig.|2Jc) (dotted line), the 
homogeneous fixed point on the middle branch is always 
unstable against the filamentary mode for L = 30, but is 
stable against Hopf bifurcations up to e < 16.43 (dotted 
line in Fig.|SJ). On the other hand, a Hopf bifurcation of 
the half-filament occurs already at e{ = 6.4, leading to 

periodic filament oscillations (Fig.QJa)) for e > el. This 
breathing filament undergoes a period doubling cascade 
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FIG. 5: Spatio-temporal spiking patterns: as in Fig.Hbut ( 
e = 13.15: periodic spiking, (b) e = 16.5: chaotic spiking. 

for e > 8.2 leading to chaos (c.f. Fig.@[b)). We note that 
the dynamic behavior is characterized by oscillations of 
varying amplitude around the unstable filamentary fixed 
point, which is typical for a chaotic breathing scenario. 

With further increase of e oscillations around the ho- 
mogeneous fixed point become important and a compe- 
tition between the two fixed points sets in. Thereby the 
breathing filament transmutes into a spiking filament as 
in Fig. Efa) for e — 13.15. Here two oscillations around 
the homogenous fixed point are followed by one oscilla- 
tion around the filamentary fixed point. This periodic 
spiking quickly becomes chaotic at further increase of e, 
yielding a chaotic spiking behaviour, which is character- 
ized by a sequence of almost homogeneous oscillations in- 
termitted by one large filamentary oscillation (Fig.^fb)). 
The resulting phase portrait is reminiscent of the Rossler 
attractor. While breathing oscillations are always close 
to the filamentary fixed point, the spiking oscillations 
are related to both the filamentary and the homogeneous 
fixed point. The complete bifurcation diagram with re- 
spect to e is depicted in Fig. We note that the homo- 
geneous limit cycle emerging from a Hopf bifurcation of 
the homogeneous fixed point (dotted line) at e = 16.43 
suppresses chaos with increasing e and finally becomes 
stable against filamentary oscillations for e > 20.2, and 
further on determines the dynamics of the system. 

Using the Benettin algorithm |32| we have calculated 
the two largest Lyapunov exponents Ai and A2 for vary- 
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FIG. 6: Bifurcation diagram of maxima and minima of the 
voltage U versus e. Thick dotted lines: homogeneous solution; 
thick dashed lines: period-one breathing orbit, thick solid 
lines: period-three spiking orbit. Parameters and load line as 
in Fig. He). 
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FIG. 7: The two largest Lyapunov exponents as a function of 
the bifurcation parameter e; load line as in Fig. |2] 

ing e (Fig. [7J. We observe that at most one Lyapunov 
exponent is positive, while the second one is zero for a 
periodic or chaotic orbit. Taking into account the third 
largest Lyapunov exponent we can then estimate an up- 
per limit for the fractal dimension of the attractor after 
Kaplan and Yorke [3^1 as Dky ~ 2.1. Although the sys- 
tem has an infinite number of degrees of freedom it is 
therefore only weakly chaotic. The global coupling is re- 
sponsible for the suppresion of extensive spatio temporal 
chaos 341 . 



IV. TIME-DELAYED FEEDBACK CONTROL 
OF CHAOTIC BREATHING AND SPIKING 

We now apply chaos control to stabilize a particular 
UPO of period t of the uncontrolled system which is em- 
bedded in a chaotic attractor. The control forces F a and 
F u in cqs. wm be designed such that they do not 

suppress chaos by generating a new periodic solution but 
only by changing the stability of a solution that already 
exists in the uncontrolled system. Hence F a and F u van- 
ish exactly on the target orbit. Such a requirement is in 
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Type of control 


F a 


F u 


diagonal control 


F\oc 


F v i 


global control with voltage feedback 


F s io 


F vf 


global control without voltage feedback 


F slo 





pure voltage control 





F v! 


local control without voltage feedback 


-Floe 






TABLE I: Overview of different control schemes with the cor- 
responding choices of F a and F u 

general met by time-delayed feedback schemes Q. For 
example, we may choose F u = F v f with a generic voltage 
feedback force defined by 



F vf (t) = u(t) - u(t - r) + RF vl {t - t). 



(4) 



Here R is a memory parameter which damps sudden 
changes in the control force by taking into account multi- 
ple delays with a decaying weight of earlier states of the 
system jjj. A control scheme of this type is called ex- 
tended time-delay autosynchronization (ETDAS). Such 
a scheme allows for noninvasive control, since an unstable 
state may be maintained with vanishingly small control 
forces. 

For the control force F a in the spatial degree of freedom 
we will alternatively use cither a local control force, 

Fio C (x, t) = a(x,t) - a(x,t - r) + RFi oc (x,t -r) . (5) 

where every spatial point is controlled independently of 
its neighboring points, or a global control force 



F s i (t) = (a)(t) - {a){t - r) + RF glo {t - r) 



(6) 



with (a) = j- adx where the same control force is ap- 
plied to every spatial point. This second choice F a = F g i Q 
may be experimentally favorable since the spatial aver- 
age is related to the total charge in the quantum well 
and does not require a spatially resolved measurement. 
In the following we will concentrate on the question how 
the coupling of the control forces to the spatial and dis- 
crete degrees of freedom influences the performance of 
the control. The UPO in question will be a flip orbit, i.e. 
its largest complex valued Floquet exponent A will obey 
ImA = ix jr. This type is practically important, as it nat- 
urally arises in period doubling scenarios, and the torsion 
of the orbit associated with ImA is a necessary ingredient 
for time-delayed feedback control to work 35] . Recently 
time-delay methods have also been extended to the sta- 
bilisation of orbits without torsion by adding a controller 
associated with an additional unstable degree of freedom 

it 

Physically, the control forces may be realized by ap- 
propriate electronic control circuits. KF U corresponds to 
an additional control voltage applied in series with the 
bias Uq, and KF a may be implemented by a spatially 
extended lateral gate electrode which influences the two- 
dimensional electron gas in the quantum well locally or 
globally. 



In principle the control performance of time-delayed 
feedback methods can be studied by linear stability anal- 
ysis of the differential-delay equations Q;© around the 
target orbit. This is difficult for a general control scheme, 
since the stability of the orbit is governed by Floquet 
eigenmodes and by the largest complex valued growth 
rates A (Floquet exponent), which are modified by the 
control force in a complicated way. Still there exists a 
particularly simple control scheme, called diagonal con- 
trol |3f||33, for which A(K, R) satisfies the exact equa- 
tion 



A + K 



1 - e- AT 



1 — Re 



-At 



= A. 



(7) 



where A is the complex Floquet exponent of the uncon- 
trolled system. Another control scheme for which an- 
alytic results are available is the Floquet mode control 
where the control force is a projection onto the unstable 
Floquet mode HHH. 

In our model we achieve diagonal control by setting 
F a = F loc and F u = F vl (cf. Table 0! , which is a straight- 
forward extension of the diagonal control for discrete sys- 
tems (cf. HH). 

We now want to control the period-one orbit (dashed 
line in Fig.^J in the chaotic breathing regime (Fig.0£b)). 
As explained in Sec. IIIII this orbit is generated by the su- 
percritical Hopf bifurcation of the stationary filament at 
e = e'j. It subsequently becomes unstable by the period 
doubling cascade with increasing e. From Fig. OJb) we 
note that, before the control is switched on (K = 0), a 
sawtooth-like structure in the (fictitious, since not yet 
applied) control force \F\ S = sup{|F v f|, |-Fioc|} appears. 
This happens as the uncontrolled system approaches and 
abandons the UPO ergodically. The slope of the increas- 
ing control force is given by A. Such a pattern is useful 
for tuning the critical parameter r empirically. After 
the diagonal control is switched on, the control force de- 
cays exponentially as \F\ S ~ | exp A(K, R)t\ (Fig. Etb)) 
to a new level which is about three orders of magnitude 
smaller than the uncontrolled level. At the same time 
the voltage signal becomes periodic (Fig.lHJa)), and the 
chaotic attractor in the phase portrait collapses to a one 
dimensional periodic orbit. 

By changing the control parameters K and R we ob- 
serve that the regime of successful control in the K- 
R plane (Fig|5{a)) exhibits a typical triangular shape, 
bounded by a flip instability (Re A = 0, Im A = tt/t) 
to its left and by a Hopf (torus) bifurcation to its lower 
right. These two boundaries may also be obtained by 
solving eq. for ReA = 0. We find that the analyti- 
cal prediction for the control boundaries are in excellent 
agreement with the numerical results (see Fig.[§fa)). 

To confirm the bifurcations at the boundaries we con- 
sider the real part of the Floquet spectrum of the orbit 
subjected to control for varying K and R = (Fig.^b)). 
Complex conjugate Floquet exponents show up as dou- 
bly degenerate pairs. The largest nontrivial exponent de- 
creases with increasing K and collides at negative values 
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FIG. 8: Diagonal control in the DBRT, where the control 
force is switched on at t — 5000. (a) Voltage u vs. time, (b) 
Supremum of the control force vs. time, (c) Phase portrait 
(global current vs. voltage) showing the chaotic breathing 
attractor and the embedded stabilized periodic orbit (thick 
solid curve). Parameters: r — —35, e = 9.1, r = 7.389, 
K = 0.137, R = 0. 
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FIG. 9: (a) Control domains in the K—R parameter plane for 
diagonal control of the unstable periodic orbit of Fig. 4(b) 
with period r = 7.389. • denotes successful control in the 
numerical simulation, ■ denotes no control, lines: analytical 
result according to eq. 0. (b) Leading real parts A of the 
Floquet spectrum for diagonal control in dependence on K 
(R = 0), • denotes analytical results from eq. J7[. 



with a branch coming from negative infinity. As a re- 
sult a complex conjugate pair develops and the real part 
increases again. The real part of the exponent finally 
crosses the zero axis giving rise to a Hopf bifurcation. 
The numerical simulations agree well with the analytical 
result. 

We now replace the local control force F a = F\ oc by the 
global control F a — F g \ without changing the voltage 
feedback (i.e. global control with voltage feedback in 
table 0. Fig. I1UI shows the corresponding control regime 
and Floquet spectrum. The control domain looks similar 
in shape as for diagonal control, although the domain for 
the global scheme is drastically reduced. The shift in the 
control boundaries is due to a different scenario in the 
Floquet spectrum. Now ReA decreases more slowly with 
increasing K than in the diagonal case, and thus the flip 
bifurcation takes place at large values of K. At the same 
time the complex conjugate pair responsible for the Hopf 
bifurcation crosses the zero axis with a larger slope and 
therefore at a smaller K value than in the diagonal case. 



It is now interesting to note that if we keep F a = F g \ a 
as before but remove the voltage feedback completely, the 
control domain is shifted to higher K values and at the 
same time is dramatically increased (Fig. Illf a)). From 
the Floquet spectrum we see that after the flip bifurca- 
tion the largest Floquet exponent does not immediately 
hybridize into a complex conjugate pair, but the Hopf 
bifurcation is caused by another complex conjugate pair 
which is not connected to the largest Floquet exponent. 
Thereby the Hopf bifurcation is suppressed and the con- 
trol regime is increased. This behavior is very similar to 
the one observed in a different reaction-diffusion model 
(modeling a heterostructure hot electron diode, HHED) 
[lq|. where it was found that additional control of the 
global variable u may gradually reduce the control regime 
to zero. 

From the practical point of view the most relevant con- 
trol scheme is the pure voltage control, i.e. F u = F v {, 
F a = 0, since the voltage variable may be conveniently 
measured and manipulated by an external electronic de- 
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FIG. 10: Same as Fig.[5]for global control with voltage feed- 
back ((b): R = 0.7) 
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FIG. 12: Same as Fig. [5]for pure voltage control ((b): R 
0.6) 



(a) 1 



0.5 



DC 



-0.5 V 




FIG. 11: Same as Fig. |5] for global control without voltage 
feedback ((b): R = 0.1). Note that the scale of the K-axis is 
changed. 



vice. The corresponding control domain and Floquet ex- 
ponents are shown in Fig. 1121 Here the control regime is 
even somewhat smaller than in the case of global control 
with voltage feedback but the shape of the control regime 
and the Floquet spectrum are qualitatively very similar. 
It is encouraging that this kind of control is at all pos- 
sible, since it opens up the opportunity to conveniently 
study chaos control in a real world device. 

We finally consider the case of local control without 
voltage feedback (Fig. I13|) . Here the control regime is 
surprisingly even larger than for diagonal control. The 
shape of the control regime is not triangular any more as 
before, but has an additional edge at low K and R val- 
ues. The reason for this edge can be explained from the 
Floquet diagram. Here a Floquet exponent from below 
collides with the largest Floquet exponent at positive real 
values before the flip bifurcation occurs. This complex 
conjugate pair then crosses the zero axes from above and 
undergoes an inverse Hopf bifurcation. For larger values 
of K another complex conjugate pair performs a second 
Hopf bifurcation seemingly unrelated to the first one. 

An overview of the different coupling schemes con- 
sidered in this section is given in table D By compar- 
ing the different control schemes, we may characterize 
the influence of the voltage control. If the voltage con- 
trol is switched on, the largest Floquet exponent rises 
quickly after the collision with an exponent coming from 
below, and participates in the Hopf bifurcation at the 
right boundary of the control regime. If on the other 
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FIG. 13: Same as Fig.[§]for local control without voltage feed 
back ((b): R = -0.55). 

hand the voltage control is switched off, the largest Flo- 
quet exponent only rises slowly after the collision, and of- 
ten the Hopf bifurcation at the right boundary is caused 
by a complex conjugate pair which is indepent of the 
largest Floquet exponent. In this case the Hopf bound- 
ary is shifted to larger values of K than for diagonal 
control. The choice of the control force F a influences 
the decrease of ReA(AT) at small K. For F a = F\ oc this 
decrease is large and the flip boundary practically coin- 
cides with the one obtained for diagonal control, while 
for F a — Fgi the flip boundary is shifted to higher K 
values since \8k ReA(AT)| is small. For F a — the slope 
| 9k Re A (AT) | is even smaller, shifting the boundary to 
even higher K values. 

The period-one orbit can be controlled continuously in 
a whole range of values of e (cf. thick dashed lines in 
Fig. HJ1 for diagonal control) . Note that since the period 
t of the UPO depends on e, t needs to be readjusted 
while sweeping e. It is then even possible to stabilize 
the period-one orbit in higher periodic windows, where 
the target orbit is obviously not part of the attractor. 
This opens up the possibility of obtaining stable self- 
sustained voltage oscillations independently of parameter 
fluctuations. The Floquet exponent A also depends on e, 
and this allows us to study the control performance as a 
function of the largest Floquet exponent of the UPO. In 
Fig. ll4f a) the control performance for global control with 
voltage feedback for fixed R is plotted. As expected, the 
regime of control is considerably smaller than the analyt- 
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FIG. 14: Control domain in the A — K plane (a) for global 
control with voltage feedback (R — 0.1), •: simulation; dot- 
ted line: analytical result for diagonal control from eq. JjJ, 
(b) pure voltage control (symbols: simulation with different 
values of R as given in the legend; lines: analytical results for 
diagonal control from eq. JjJ). 



ical predictions for diagonal control. Fig. ll4f b) shows the 
control performance for pure voltage feedback and differ- 
ent values for R. Although the control regimes are much 
smaller than predicted for diagonal control, the trend in 
the shift of the control regime for increasing R is quali- 
tatively similar. 

Surprisingly we find a finite K value for the flip bound- 
ary as At — > (Fig. IHT b)). This is in marked contrast 
to the analytical results for diagonal control, where the 
flip boundary extrapolates to zero for At — > 0. To un- 
derstand why this is the case, let us consider a situation 
where we apply pure voltage control to a stable orbit. 
From the Floquet spectrum in Fig. 1151 we see that with 
increasing control force the stable orbit is destabilized by 
a flip bifurcation, then becomes stable again by a second 
flip bifurcation, before it finally undergoes the usual Hopf 
bifurcation at high K values. Such a reentrance scenario 
does not happen for diagonal control. 

So far we have stabilized the period-one breathing or- 
bit. We have investigated whether the obtained results 
are specific for a given UPO, or whether similar results 
can be observed for different orbits such as the period- 
three spiking orbit (cf. thick solid lines in Figs. l5ltj|l which 
stands for the second basic dynamical pattern (spiking) 
in the DBRT system. Fig. ^| shows the control regimes 
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FIG. 15: Leading real parts A of the Floquet spectrum for 
pure voltage control of a stable period-one breathing orbit in 
dependence on K (R — 0.6, e = 8.0) 
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FIG. 16: Control domains for spiking orbit. Lines: analytical 
results for diagonal control, (a) diagonal control, (b) global 
control without voltage feedback, (c) local control without 
voltage feedback (Parameters: e = 13.44, r = 20.15) 

for different control schemes. Again we note that di- 
agonal control is very well predictable by the analytic 
formula (7J (Fig. IlKf a)). whereas omitting the voltage 
control shifts the Hopf border to large K values, and 
replacing local with global control shifts the flip bound- 
ary to larger K values. For pure voltage control and 
global control with voltage feedback no domain of control 
was found. Nevertheless features of the different control 
schemes are similar for both breathing and spiking orbits. 

V. CONCLUSIONS 

We have applied different schemes of chaos control 
by time-delay autosynchronization to a globally coupled 



reaction-diffusion model describing charge transport in 
a semiconductor nanostructure, viz. a double-barrier 
resonant tunneling diode (DBRT). The spatio-temporal 
dynamics and bifurcation scenarios of the uncontrolled 
DBRT in an external circuit display a variaty of complex 
patterns, including breathing and spiking oscillations ei- 
ther of which can be periodic or chaotic. 

We have shown that by using time-delayed feedback 
methods it is possible to stabilize unstable breathing or 
spiking patterns. Different control scheme were com- 
pared with respect to efficiency, and quantitative compar- 
isons of their respective control domains in the R — K pa- 
rameter space were given and interpreted in terms of the 
Floquet spectra. Those turned out to be helpful to gain 
insight into the control mechanisms and understand why 
the control performance may, e.g., be improved by omit- 
ting control of the voltage variable. As a control scheme 
which is particularly simple to realize experimentally, al- 
beit with a small control regime which requires careful 
adjustment of the control parameters R and K, we have 
singled out pure voltage control. Our findings may be 
important for obtaining stable self-sustained voltage os- 
cillations in resonant tunneling diodes independently of 
parameter fluctuations. 

This work was supported by DFG in the framework of 
Sfb 555. 



APPENDIX A 

In this Appendix we provide analytical expressions for 
the functions D(a) and f{a,u) appearing in eq. and 
relate the dimensionlcss quantities used throughout this 
paper to their respective dimensional physical quantities 
marked by a tilda. 

The voltage drop across the DBRT u(t) and the 
external bias voltage Uq are related to their dimen- 
sionless counterparts by u ~ e-u/(fcsT) and Uq = 
ello/ '(ksT), respectively, with the temperature T, the 
electron charge e < 0, and Boltzmann's constant ks- 
The two-dimensional electron density is rescaled via 
a = airh 2 / '{rnksT), where m is the effective mass of 
the electron in the well. Time and sp ace coordinates 
are rescaled by t — tY^lh and x = x / ^h^iksT / (— eT^) 
where Tl = Tr is the transition rate for electrons from 
the emitter to the quantum well (and from the quantum 
well to the collector) and /i is the electron mobility. 

This yields a rescaling of current and resistance as j = 
jttH 3 / (emkBTT l) and r = re 2 mAT l / (ttH 3 ) , respectively, 
where A is the cross sectional area of the device. Note 
that the effective resistance in general contains two terms 
f = R — k, where R represents an external resistor and 
k arises from an active circuit with a voltage gain U c — 
kjA [23. In the case k > R this leads to r < with 
particularly interesting chaotic scenarios, as shown in the 
main text. The time-scale ratio is given by e = RCTl/H, 
where C is the total capacitance of the circuit. 

The effective diffusion coefficient D(a) may be calcu- 
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lated by considering a generalized form of Einstein's re- 
lation which covers the full range from Fermi-degenerate 
to nondegenerate conditions and a drift term resulting 
from the change in the local potential due to Poisson's 
equation |2{j: 



D(a) = a — 

JB 



1 



1 — exp [—a] 



(Al) 



where tb = (47reeofi- 2 )/(e 2 m) is the effective Bohr radius 
in the semiconductor material, e and eo are the relative 
and absolute permittivity of the material, respectively, 
and d is the effective thickness of the double-barrier struc- 
ture. 

The function f(a,u) is obtained from a microscopic 
consideration of the tunneling currents J ew (a,u) and 
Jwc(a>) from the emitter to the quantum well and from 



the well to the collector, respectively [2E 



f(a,u) = 
In [ 1 + exp (r\ e 



• arctan 



-(xo 



7 



u 
2 



X 



u d ■ 

77 a ) — a 

2 r B ' 1 



(A2) 



Here 7 and Xo describe the broadening and the energy 
level of the electron states in the quantum well, and r) e is 
the dimensionless Fermi level in the emitter, all in units 
of UbT . In this article we choose 7 = 6, djrs = 2, r\ e = 
28 and £0 = 114. Typical physical values correspond 
to units of time, space, voltage, electron density, and 
current density of the order of 3.3 ps, 100 nm, 0.35 mV, 
10 10 cm -2 , 500 Acm~ 2 , respectively, for d = 20 nm, — 
T R = 0.2 meV, k B T = 0.32 meV (T = 4 K ). 
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